UNCLASSIFIED 


AD  NUMBER 


AD479489 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution:  Further  dissemination  only 
as  directed  by  Office  of  Naval  Research, 
Arlington  VA  22217-0000;  Sep  1965  or 
higher  DoD  authority. 


AUTHORITY 


ONR  ltr  27  Jul  1971 


THIS  PAGE  IS  UNCLASSIFIED 


¥7f  Vff 


! 


HYDRODYNAMICS  LABORATORY  I 


KARMAN  LABORATORY  OF  FLUID  MECHANICS  AND  JET  PROPULSION 


CALIFORNIA  INSTITUTE  OF  TECHNOLOGY 


PASADENA,  CALIFORNIA 


Hydrodynamics  Laboratory 

Karman  Laboratory  of  Fluid  Mechanics  and  Jet  Propulsion 
California  Institute  of  Technology 
Pasadena,  California 


EVALUATION  OF  PRESSURE  DISTRIBUTION 
ON  A  CAVITATING  HYDROFOIL  WITH  FLAP 

by 

Zora  L.  Harrison 
and 

Duen-pao  Wang 


This  work  was  carried  out  under  the 
Bureau  of  Ships  General  Hydrodynamics  Research  Program, 
Administered  by  the  David  Taylor  Model  Basin, 

Office  of  Naval  Research  Contract  Nonr  -220(52). 


Reproduction  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government 


Report  No.  E-133.  1 


September  1965 


ABSTRACT 


A  general  method  is  established  to  calculate  the  pressure  distri 
bution  and  the  moment  of  force  for  a  two-dimensional,  supercavitating 
hydrofoil  with  a  flap.  The  wake  flow  model  is  adopted  to  describe  the 
configuration  of  the  flow  field.  Some  numerical  results  for  a  super¬ 
cavitating  flat  plate  with  a  flap  are  compared  with  the  corresponding 
experimental  data. 


NOMENCLATURE 


A 

a  ,  a 
1  2 

AR 

c  ,  c 
1  2 


a  scale  factor 

distances  from  the  hinge  point,  P,  to  the  respective  centers 
of  pressure  of  the  main  body  and  the  flap 

aspect  ratio 

distances  along  the  chord,  measured  from  the  hinge  point 
total  drag  coefficient 

total  normal  force  coefficients  acting  at  the  centers  of  pressure 
total  lift  coefficient 

moment  coefficients  about  the  hinge  point 

moment  coefficient  of  the  flap  about  the  hinge  point,  CMtt=  CM 
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total  moment  coefficient  about  the  leading  edge 
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force  components  in  y-direction 

complex  potential 
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chord  lengths:  =1,  f/c  =  i  /  (i^  +  i  ) 

function  in  the  integrand  of  V 

& 

distance  parameters  in  the  calculation  of  and  (3^ 

(note  the  subscripts  l  and  2  do  not  refer  to  the 
main  body  and  the  flap  for  these  parameters) 

the  prescribed  value  of  f/c 

chord  distance  measured  from  the  leading  edge 

a  mapping  plane  of  the  flow  region 

image  point  of  z  -•  oo  in  the  t-plane:  t  =  Ve~*^ 

image  of  the  hinge  point  in  the  t-plane 

-2 

free  stream  velocity:  U  =  1  +  o 
magnitude  of  t 
complex  velocity 

coordinate  in  z -plane  and  transformation  variable 
coordinate  in  z-plane 
physical  plane 
angle  of  attack 

polar  angle  of  t  in  the  t-plane 
flap  angle 

angles  between  the  tangents  to  the  foil  profile  and  the  chords 

cavitation  number 

integration  variable  in  the  t-plane 
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Introduction 


In  practical  applications  of  supercavitating  hydrofoils  the  adoption 
of  flaps  or  other  load  modulators  as  motion  control  devices  is  of  funda¬ 
mental  importance.  An  accurate  determination  of  the  moment  and  hinge- 
moment  acting  on  supercavitating  hydrofoils  with  deflected  flaps  is  essen¬ 
tial  to  the  structure  design  on  one  hand,  and  to  the  development  of  success¬ 
ful  control  and  maneuvability  on  the  other.  Furthermore,  a  detailed  study 
of  the  pressure  distribution  is  also  valuable  for  establishing  the  general 
validity  of  any  theoretical  model  and  also  for  elucidating  the  role  of  effec¬ 
tiveness  of  the  flap. 

Recently,  the  supercavitating  flap  problem  has  been  investigated 
using  the  linearized  theory  for  the  case  of  zero  cavitation  number,  corres¬ 
ponding  to  an  infinitely  long  cavity.  The  case  of  infinite  flow  extent  was 
first  treated  by  Tulin  and  Burkart  [l]  .  This  theory  has  been  extended  by 
Johnson  [2]  who  incorporated  the  classical  work  of  Kirchhoff  and  Rayleigh, 
and  that  of  Green,  the  latter  including  the  effect  of  the  free  water  surface, 
but  neglecting  the  gravity  effect.  In  Johnson's  tneory,  an  estimate  of  the 
effect  of  finite  aspect  ratio  is  also  evaluated.  A  linearized  theory  was 
subsequently  developed  by  Auslaender  [3]  for  the  case  of  zero  cavitation 
number  and  at  a  finite  depth  of  submergence.  Using  a  nonlinear  theory, 

Lin  [4]  treated  the  flap  problem  for  the  flat  plate  hydrofoil  at  zero  cavita¬ 
tion  number.  In  an  approximate  sense,  the  limiting  condition  of  zero 
cavitation  number  corresponds  to  the  shallow  submergence  case  when  the 


1  Numbers  in  brackets  designate  references  at  end  of  paper. 


cavity  become?  completely  open  to  the  atmosphere.  However,  a  fine 
difference  in  this  correspondence  arises  from  the  fact  that  the  dynamic 
effect  due  to  the  upper  half  flow  field  is  actually  modified  in  the  pre¬ 
sence  of  a  free  water  surface  at  shallow  submergences.  In  accordance 
with  this  analogous  configuration,  some  model  investigations  with  flaps 
have  been  carried  out  by  Johnson  [2]  ,  Brown  [5]  ,  Conolly  [6]  and  Wetzel 
and  Maxwell  [7]  .  In  order  to  establish  the  fully  cavitating  configuration 
at  small  incidences,  the  technique  of  base  ventilation  has  been  employed 
in  some  cases  of  these  experiments. 

The  objective  of  the  present  investigation  is  to  establish  a  general 
and  systematic  method  for  calculating  accurately  the  pressure  distribu¬ 
tion,  moment  of  force,  and  hinge -moment  acting  on  a  supercavitating 
hydrofoil  of  arbitrary  profile,  especially  when  the  angle  of  attack,  flap 
deflection  and  the  cavitation  number  are  moderate  or  large.  This  calcu¬ 
lation  is  based  on  the  nonlinear  theory  developed  recently  by  Wu  and 
Wang  [8]  .  To  meet  the  general  operating  conditions  in  practice,  this 
theory  can  be  further  extended  to  include  the  effects  of  free  water  sur¬ 
face  and  gravity,  and  finally  the  effect  of  finite  aspect  ratio  and  three-dim¬ 
ensional  features.  The  scope  of  the  present  work,  however,  will  be  limited 
to  the  case  of  two-dimensional  flow  of  infinite  extent,  with  emphasis  on  the 
pressure  distribution  over  a  flapped  hydrofoil.  The  formulas  for  the  force 
moment  and  hinge -moment  are  also  derived  here.  Since  the  information 
about  pressure  distribution  in  cavity  flows  is  relatively  scarce  in  the  ex¬ 
isting  literature,  it  is  hoped  that  the  resul'i  of  this  study  may  provide  use¬ 
ful  data  for  future  applications  as  well  as  for  further  development  of  this 
general  field. 


The  numerical  computation  has  been  performed  with  an  IBM  7094 
Computer,  for  which  the  method  and  computing  program  are  explained 
in  detail  in  this  report.  The  final  results  for  the  case  of  flat  plate  hydro¬ 
foil  with  a  flap  are  presented  in  figures  'or  different  angles  of  attack, 
flap  angles,  flap-chord  ratio  and  the  cavitation  number.  The  calculation 
for  hydrofoils  of  arbitrary  profile  is  also  discussed.  The  computation 
program  for  this  case  fo.llow3  an  approximate  numerical  scheme  developed 
by  Wu  and  Wang  [9]  .  This  approximate  method  greatly  shortens  the  com¬ 
putation  of  the  final  integral  equations  in  the  theory,  while  still  retaining 
a  high  degree  of  accuracy  of  the  result. 

A  parallel  experimental  program  has  been  accomplished  by  Meijer, 
the  results  of  these  experiments  will  be  presented  in  a  separate  report  by 
Meijer  [lO],  It  can  be  stated  here  that  the  general  agreement  between  the 
theory  and  experiments  is  very  satisfactory,  and  the  experimental  obser¬ 
vation  has  also  revealed  some  interesting  points,  such  as  the  viscous 
effect  near  the  flap  hinge.  With  the  theory  now  well  established  in  the 
typical  cases  so  far  investigated,  the  present  method  can  be  applied  with 
confidence  for  the  g'-r>eral  case  of  arbitrary  profile. 

General  Formulation 

We  consider  a  steady,  irrotational  plane  flow  of  an  incompressible 
fluid  past  an  arbitrary  profile  AB  in  such  a  way  that  the  flow  is  separated 
from  the  leading  and  trailing  edges  A  and  B,  forming  a  fully  developed 
cavity  behind  the  foil  AB.  The  foil  AB  consists  of  a  main  foil  AP  and  a 
flap  PB,  as  shown  in  Fig.  1(a).  A  set  of  Cartesian  coordinates,  x  and 
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y,  is  chosen  with  the  origin  at  the  leading  edge  A  and  the  x-axis  lying 
along  the  chord  AP.  The  incoming  free  stream  has  a  speed  U  and  is 
at  an  angle  of  attack  a  to  the  x-axis.  The  flap  angle,  t  ir,  is  the  angle 
measured  from  the  x-axis  to  the  chord  of  the  flap  PB,  positive  in  the 
clockwise  sense. 

In  the  following  analysis  the  wake  flow  model,  as  formulated  in 
[8]  ,  will  be  used  to  describe  the  flow.  For  the  sake  of  convenience,  a 
general  description  of  the  wake  model  will  be  given  briefly  in  the  follow¬ 
ing:  The  part  AC  and  BC1  of  the  free  streamlines,  as  shown  in  Fig.  1(a), 
form  the  boundary  of  a  near -wake  of  constant  pressure  pc,  which 
is  less  than  the  free  stream  pressure  From  this  part  onward  the 

pressure  varies  continually  and  monotonically  from  pc  to  p^,  along  the 
far -wake  boundary,  CX  and  C'l.  It  is  further  assumed  that  the  values 
of  the  complex  potential,  f ,  and  the  complex  velocity,  w ,  at  the  point 
C  are  equal  to  those  at  the  point  C'  respectively.  With  the  introduction 
of  the  so-called  "hodograph-slit  condition"  for  the  streamlines  Cl  and 
C'l,  we  assume  that  the  entire  flow  region  in  the  z  -plane  can  be  mapped 
onto  the  interior  of  the  lower -half  unit  circle  in  a  1-plane,  as  shown  in 
Fig,  l{b),  where  T  is  the  image  of  the  point  P,  and  t  that  of  the 
point  I.  On  the  t-plane  the  complex  potential,  f,  and  the  complex  veloc¬ 
ity,  w,  with  the  speed  on  the  cavity  wall  normalized  to  unity,  can  be 
written  respectively  as 
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where  6  (z(t))  is  the  angle  measured  from  the  chord  AP  (or,  chord  PB), 
counter-clockwise  as  positive,  to  the  tangent  at  any  point  z  on  the  sur¬ 
face  profile  curve  AP  (or,  curve  PB),  as  shown  in  Fig.  1(a).  At  z  =  ■», 


w  =  tie  and  t  =  t  ,  hence 
o 


For  any  point  on  the  foil,  corresponding  to  t  *  t,  t  being  real,  -  1<t<  1, 
the  distance  along  AP  and  PB  between  the  leading  edge  A  and  the  pro¬ 
jection  of  this  point  to  AP  or  PE  is  given  by: 

c(t)  =  2a£ f  g(t;tQ)  |exp[if£ 


where  *  above  the  integral  sign  signifies  that  the  Cauchy  principal  value 
of  the  integral  is  taken  and 
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[(t-to)  <t  - to)  (t-t/1)]* 


The  condition  that  the  points  G,  C'  lie  on  the  lower  unit  circle  in  the 
t-plane,  or  -ir<argtc<0,  may  be  expressed  as 

-1  <Re  i(tQ,  —  )<  1  . 
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This  condition  is  equivalent  to  stating  that  the  constant  pressure  region 
extends  beyond  the  trailing  edge  B. 

By  writing  t  =  Ve~^,  Eq.  (4)  becomes 

(t;t  }  -  V4(l  -t*)[l+<*-(l+V~2)(tV  cos  p)]  (5 

°  [(t2  +  V2-  2tV  cos  p  Mt2  V2  -  2 tV  cos  p  +  1 )]  2 

Similarly,  substituting  t  =  Ve-i^  into  Eq.  (2),  expanding  and  separating 
the  real  and  imaginary  parts,  we  obtain  the  parametric  relations  for  t 
as  follows 
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Q(t)  = 


tv2  -  V  cos 


t2  V2-  2tV  cos  p  +  1 


t -  V  cos  (3 
t2  +  V2  -  2tV  cos  p 


H{t)  =  - - -  +  - - -  . 

t2  V2 -2tV  cos  p  +  1  t*  +  V2  -2tV  cos  P 

L.et  J{t)  denote  the  Cauchy  integral  in  Eq.  (3).  Then,  noting  that 


f 1  dr 
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the  singularity  o£  the  integrand  in  J (t)  can  be  removed,  giving 


J(t)  ,  f  *  Wd T  =  C*  0(T)  -  fl  (t) 
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If  we  let  the  chord  of  AP  be  i  and  the  chord  of  PB  be  i.  then 

1  2 

=  2A  J  I  ~T~Tt  '  1  8(*  i  t0 )  J^tjJ  cos6(t)}dt 


=  2A  C  j  1  g{t;  tQ)|exp£^— J^(t) J  cos  6(t)|  dt 


where  i  ,  l  ,  6  and  9  are  shown  in  Fig.  lc.  For  given  U,  a,  9  lx), 
1  2  1  2 

and  1^,  Eqs.  (2),  (3)  and  (17)  provide  a  set  of  functional  equations 


f  The  subscripts  1  and  2  in  all  following  equations  refer  to  the  main  foil 
and  the  flap  respectively. 
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for  unknown  t  ,  c(t),  A  and  T. 

To  facilitate  the  numerical  computation  of  the  above  integrals  it 
is  convenient  to  introduce  the  following  changes  of  variables 

iw.ii,  ,  W,JE±*L;  (18) 

1  1  -  Tx2  2  1  +  Tx2 


so  that  t  varies  from  T  to  -1,  t  from  T  to  +1,  as  x  increases 
1  2 

from  0  to  1.  Then  Eq.  (17)  may  be  written  as 

ik  =  4A(1  -T2)^  (19) 

where,  with  k  =  1,2, 
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2  2  2  2 

6  (x)  and  0  (x)  in  Eqs.  (20)  and  (21)  are  the  angles  between  the  tangents 
1  2 

and  the  chords  of  AP  and  PB  at  c^  (x)  and  c^  (x)  respectively,  as 
shown  in  Fig.  1(c). 

If  we  now  normalise  the  total  chord  (1^  +  )  to  be  unity,  we  have 
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and  hence 
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The  ratio  of  the  flap  chord  f  to  the  total  chord  c  is 


(23) 


f/c  =  i4-  -  rrr 

12  12 


(24) 


Making  the  same  changes  of  variable,  Eqs.  (9)  and  (11)  become 
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where  Q  and  H  are  given  by  Eqs.  (14)  and  (15)  respectively. 
The  pressure  coefficient  at  any  point  is 
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P  ipU2 


(27) 


where  p  is  the  pressure  at  the  wetted  surface  of  the  plate  and  pc  the 
cons  tent  pressure  in  the  cavity.  From  the  Bernoulli  equation  we  obtain 
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Cp  =  (I+ff)(l  -w  w  )  , 


(28b) 


where  cr  is  the  wake  under -pressure  coefficient,  or  the  cavitation  number 

p«r  pc  — 

for  cavity  flows,  given  by  a  =  — -  ,  and  w  is  the  complex  conjugate 


ipu* 


of  w. 


The  complex  force  coefficient  for  the  main  body,  Cx  +  iC  ,  can 
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be  written  as 
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From  Eqs.  (I),  (4),  (22)  and  the  transformation  (18),  we  can  write  the 
above  equation  as 


C  +  iC  =  i(l  +  cr) 
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Similarly,  the  complex  force  coefficient  for  the  flap,  Cx  +  iC  ,  be- 
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The  lift  and  drag  coefficients,  and  ,  are  given  by 

CD+  iCL*  t(Cx  +  iCv  >  +  (cx  +  iC  )i  e_ia  •  (31) 

1  yi  2  * 2 

Special  Case:  Flat  Plate 

As  a  special  case,  we  now  consider  a  flat  plate  with  a  flap  as 
shown  in  Fig.  2.  This  configuration  is  particularly  simple  since  6=0 
everywhere  along  the  chord.  By  putting  6  =  0  in  our  formulae,  V 
given  by  Eq.  (9)  becomes  unity  and  jJ  given  by  Eq.  (1 1 )  becomes  zero 

3L 


so  that 

> 

ii 

> 

(32) 

and 

-G> 

II 

(33) 

where  and  are  given  by  Eqs.  (8)  and  (10)  respectively.  Similarly 
w,  3y  1^,  ^  and  c^  are  given  by  Eqs.  (lb),  (20)  and  (23)  with  the 

factor  inside  the  curly  brackets  equal  to  unity. 

From  the  above  and  Eq.  (28b)  we  obtain  the  pressure  coefficient 

Cp  =  (l  +  a)[l -t2  |^-t  I”]  •  <34) 

In  terms  of  the  variable  x  defined  by  Eq,  (18),  Eq.  (34)  may  be  written 
as 

Cp^(x)  a  (1  +<r)  [l  -tj^2  (x)x4*  ]  ,  (k  a  1,2).  (35) 

The  moment  coefficients  about  the  hinge  point  P  due  to  the  pressure 
over  the  main  foil,  AP,  and  of  the  flap,  PB,  as  shown  in  Fig.  2  ,  are 

f  T 

C..  =  \  C  (t)  c  (t)  dc 

Mi  Pi  1  1 
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and 

cm2=  rv  c2(t) 

respectively.  From  Eq.  (23)  we  get 

l  T1  x(!-20  gOUxbt  )fx  (l-2t)  gtUytJt  ) 

+  I^)2  4t>  pk  [l  +  (-r Tx2]2  V4  [l +  (-TTy2]4  V4 

(k  =  1,2)  (37) 

where  1^  and  are  given  by  Eq.  (20)  with  the  factor  inside  the  curly 
brackets  equal  to  unity  and  g(t;  tQ)  is  given  by  Eq.  (4). 

It  is  also  desirable  to  know  the  location  of  the  centers  of  pressure 
of  both  the  main  body  and  the  flap  as  well  as  the  total  forces  at  those 
points.  Let  and  a^  be  the  distances  from  the  hinge  point,  P,  to  the 
center  of  pressure  on  the  main  foil  and  to  that  on  the  flap  respectively, 

and  let  and  C^,  be  the  force  coefficients  acting  at  these  points. 

i  z 

We  may  now  write 

CFi=I,lCpdC‘  °V  CpdC  ’ 

with  the  help  of  Eq.  (23)  the  final  results  are 


x<*  -  *) 

[l  +  <-)kTx2F 


8{tk(x);to}  , 
- dx, 

V4 


(k  »  1,2), 


(38) 


*>■  ■  cvc^k  ■ 


(k  =  1,2). 


(39) 
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Knowing  these  quantities  we  can  transfer  the  total  moment  to  the  leading 
edge,  giving 

CM  "  CF(Vax)+  CF  co»*ir  +  a2),  (40) 

LE  i  2 

where  C.,  is  considered  to  be  positive  in  the  nose  down  sense.  The 
MT,1T 

hinge  moment  is  defined  as  the  moment  of  the  flap  about  the  hinge  point, 
hence 


M 


■  C. 


H 


M  ‘ 
2 


(41) 


The  pressure  distribution  is  given  in  the  parametric  form  s(x), 
where  s(x)  denotes  the  length  al 
from  the  leading  edge  such  that  0  <  s  <  1 


C’  ■(*),  where  s(x)  denotes  the  length  along  the  foil  surface  measured 


s(x)  =  1  - - c  (x)  on  AP 

c  1 


=  1 - +  c  (x)  on  PB 

C  z' 


(42) 


and  C1  is  the  pressure  coefficient  which  is  defined  as  unity  at  the  stagna¬ 
tion  point  D  and  is  equal  to  negative  0  on  the  cavity  side.  Then  we  have 

C*p  (x)  =  (1  +  0)  [l  -t^xjx*]  -  a,  (k  =  1,2).  (43) 

XV 

The  lift  and  drag  coefficients  are  given  by  (refering  to  Eq.  (78) 
of  Ref.  (8)) 


C 


L. 


ttA 

U2 


(g-tU)ccp  f  ,  , 

(v"2+  V2-  2  cos  2J3)  | 


-it  m2  »  T  <  1  .  m2  .  .  ..  O  m« r  f  i  i  -  -  '5(3  1 

Mi  -  i  mit  i  yv  CoS  p-j.v\iTv  f  v;QS  t*pj  \ 

(T2  +  V2  -  2tV  cos  p)(T2V2+  1  -  2tV  cosp) 


(44) 
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c  _  trA  (U'1-  U)/(l  -Vz)  L  ,  va,  *(1  -  T*)[(l  +  T2)V2(1  +  V2)  -  2TVcoaP(H-  V4)] 
D  U2  (V‘2+ Vz-2cos2p)  \  <T2+ V2-2TVcosp)(T2V2+ 1  -  2TVcosP) 

(45) 

where  A  is  given  by  Eq.  (22). 

Numerical  Procedures: 
a)  Flat  plate  with  a  flap. 


A  double  iteration  procedure  is  required  to  evaluate  the  force  co¬ 
efficients  for  prescribed  U,  a,  e  and  £/c.  We  first  solve  the  inverse 
problem:  choose  U,  a,  t  and  T  as  the  independent  parameters,  solve 
for  V  and  P  (hence  tQ)  by  iterating  the  parametric  relations  in  Eqs.(32) 
and  (33).  We  then  obtain  the  corresponding  flap-to-chord  ratio,  f/c 
a  +  1^),  where  1^  and  1^  are  given  by  Eq.  (20)  with  the  factors  in¬ 
side  the  curly  brackets  equal  to  unity.  A  type  of  iteration  similar  to 
Newton's  method  is  then  used  on  T  to  find  V  and  p  corresponding  to 
the  desired  value  of  f/c. 

Once  V  and  P  are  established,  the  moment  coefficients,  CM 

:le 


and 


are  obtained  f 


from  Eqa.  (42)  and  (43)  ai 


"qs.  (37)-(41),  the  pressure  distribution 
!t  and  drag  coefficients  from  Eqs.  (44) 


and  (45). 


The  details  of  the  cau  dons  outlined  above  are  described  in 
Appendix  A*. 


w 


The  program  listing  is  available  upon  request. 


b)  Arbitrary  profile. 


The  iteration  scheme  for  the  arbitrary  profile  is  a  combination 
of  the  above  method  and  the  approximate  scheme  of  Wu  and  Wang  [  9] 
and  is  described  in  detail  in  Appendix  B  . 

Results: 


The  numerical  calculations  outlined  above  were  carried  out  for 
three  models  of  flapped  flat  plate  hydrofoils  with  flap  to  chord  ratios 
of  9.  2,  0.  4  and  0.  6.  For  each  model  the  flap  angle,  c  xr,  was  varied 
from  10°  to  60°  and  the  angle  of  attack,  a,  from  5°  to  60°  for  cavitation 

^C>J{ 

numbers  between  0  and  1.0 

Figures  3  through  5  show  the  total  moment  coefficient  about  the 
leading  edge,  ,  and  the  moment  of  the  flap  about  the  hingepoint. 


LE 


C*,  ,  as  functions  of  a  for  a  =  5°  through  60°.  For  each  model  the 


M. 


H 


flap  angles  10°,  20°,  40°  and  60°  are  shown.  CM  as  a  function  of  eir. 


mt,f 

is  given  in  Fig.  6  for  all  three  models  at  various  cavitation  numbers 
for  a  =  10°,  20°  and  30°.  Figure  7  shows  the  lift  coefficient.  C^, 
versus  «  tv  for  the  same  conditions.  The  polar  plots  of  C^/Cq  versus 
CT  for  os  0  are  shown  in  Fig.  8,  in  which  the  solid  lines  are  lines  of 
constant  a  and  the  dashed  lines  are  those  of  constant  <n.  Since,  within 
the  range  of  the  cavitation  numbers  concerned,  C^/C^  is  almost  inde¬ 
pendent  of  o,  no  corresponding  curves  are  needed  for  other  values  of  o. 
Some  typical  pressure  distribution  curves  for  the  model  with 


The  program  listing  is  available  upon  request. 
Tabulated  results  are  available  upon  request. 


-■ _ so.-. 
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£/c  =  0.2  at  a  =  10°  and  «  v  =  0°,  20°,  40°  and  60°  are  shown  in  Pig.  9. 
The  experimental  points  are  taken  from  the  work  by  Meijer  [  1 0] .  In 
Meijer's  experiments,  different  sized  models  were  tested  in  order  to 
make  wall  effect  corrections.  In  Fig.  9  the  theoretical  results  agree 
very  well  with  the  experimental  ones,  except  around  the  hinge  region. 

The  discrepancy  there  may  be  attributed  to  viscous  effect.  Curves  were 
faired  through  Meijer's  experimental  points,  given  in  Fig.  9.  and  then 
integrated  to  obtain  and  C^.  The  results  are  plotted  along  with  the 
corresponding  theoretical  lift  and  drag  curves  in  Fig.  10.  As  would  be 
expected  from  the  pressure  distribution  correlation  in  Fig.  9,  the  agree¬ 
ment  in  CL  and  i  very  good. 

Figure  11  illustrates  the  agreement  between  Parkin's  experiment 

[ll]  and  the  theory  of  C. .  as  a  function  of  o  for  a  flat  plate  at  angles 

™LE 

of  attack  from  10°  to  30°.  The  experimental  points  of  Parkin's  experi¬ 
ment  have  not  been  corrected  for  wall  effects. 

Besides  the  above  comparisons  between  two-dimensional  results, 
it  is  also  interesting  to  compare  the  present  two-dimensional  theoretical 
results  with  experiments  having  three-dimensional  effect.  Wetzel  and 
Maxwell  [?]  performed  experiments  with  ventilated  finite  aspect  ratio 
models  at  different  submergences  in  a  free  jet.  Some  of  these  results 
are  shown  in  Figs.  12  and  13.  In  general,  the  agreement  between  the  two- 
dimensional  theory  and  three-dimensional  flow  is  not  as  good  as  the  agree¬ 
ment  with  the  two-dimensional  experiments  of  Parkin  and  Meijer,  but 
as  aspect  ratio  and  submergence  are  increased  the  agreement  tends  to 
improve. 


IV- 


1? 
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APPENDIX  A 


The  double  iteration  procedure  outlined  in  the  text  for  the  flat 
plate  problem  presents  no  problem  when  done  on  a  high  speed  computer. 

All  calculations  presented  here  were  done  on  the  IBM  7094  at  Booth 
Computing  Center,  California  Institute  of  Technology. 

Choose  U,  g,  (  and  let  the  first  approximation  of  T  be  0.5, 
then  calculate  V  and  P  by  iteration  of  Eqs.  (32)  and  (33).  Since  V 
and  p  are  rather  close  to  U  and  a  in  value,  let  V^=  U  and  P^  =  a 
be  the  first  approximations,  then 

V(1)  =  Vf(T,V(0},  P(0)) 

and 

P(1)  =  Pf(T,V(0),  P(0)> 

.  (n  -  U 

are  the  next  approximations.  Continue  this  procedure  until  |  V 
-  |  and  |  p^n”  ^  -  p^  j  are  less  than  the  allowable  error.  This 

iteration  converges  very  well ;  on  the  average  about  6  iterations  will 
produce  an  accuracy  such  that  the  error  is  less  than  0.  00001. 

Once  V  and  p  are  determined  the  corresponding  flap  to  chord 
ratio,  f/c,  can  be  calculated  from  Eq.  (24),  where  1^,  and  1^  are 
given  by  Eq.  (20)  with  the  curly  brackets  equal  to  one.  Let  the  desired 
f/c  be  R  and  the  first  approximations  to  T  and  f/c  be  and 

respectively.  Since  f/c  decreases  as  T  increases,  the  second  approxima¬ 
tion  to  T  is  arbitrarily  set  to  be  either 


T^1  *  =  0.  1  ;  R  -  r(°*>  0 


20 


or 

.  o  ^  *  H.  o 

The  corresponding  R^  is  calculated  from  the  first  iteration  process 
and  now  we  can  apply  a  modified  Newton's  method. 


and 

R<2>=|.  <T<2>>. 

Continue  these  successive  approximations  of  T  according  to 


until  jR-R^n-1M  is  less  than  the  allowable  error. 

It  should  be  pointed  out  here  that  only  positive  values  of  T  are 
of  interest  for  the  physical  application  considered  here.  If  X  becomes 
negative  then  the  stagnation  point  D  has  moved  to  the  flap  side  of  the 
hinge  point.  This  situation  occurs  when  f/c**l  and  a  and  c  become 
large. 

Once  T,  V,  and  p  are  determined  the  pressure  distribution, 
lift,  drag  and  moment  coefficients  are  readily  calculated  as  outlined 
in  the  text.  The  numerical  integration  required  was  programed  using 
Simpson's  rule  on  increments  of  x  ■  0.  005  and  y  =  0.  0001.  This  in¬ 
crement  size  was  sufficiently  email  to  produce  an  accuracy  c-f  at  least 
±  0.  0001  for  all  cases  that  were  calculated.  A  good  check  of  the  accuracy 
is  to  print  out  s^x  =  1)  and  s^(x  =  1)  as  these  values  should  be  G  and 
1.  0  respectively. 
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APPENDIX  B 

Th«*  iteration  procedure  required  for  the  arbitrary  profile  is 
s*im«>v/b»  f  lonj tt'r  and  more  complex  than  the  flat  plate  problem  be - 
*. »»!«*•  of  lb*'  added  elope  function,  0.  This  additional  parameter  gives 
riao  |o  an  inleg ration  in  the  V  and  p  calculations;  furthermore,  the 
rliortl  <ra lr;ola.  jt»n  liw:mnpR  a  rloublc  integral.  We  have  used  a  modifica- 
li«tii  of  Ihe  approximate  numerical  scheme  originated  by  Wu  and  Wang  [lO]  . 
This  modification  is  necessary  because  of  the  additional  requirement  of 
fixed  f/c  so  that  T  cannot  be  prescribed. as  it  was  in  that  reference. 

We  first  solve  lire  flat  plate  problem  for  the  prescribed  a,  a,  e 
and  f/c.  This  solution  gives  us  a  good  first  approximation  for  the  un¬ 
known  parameters  V.  P  and  T.  Using  these  initial  values  we  now  iterate 
l",qs.  (h).  (7)  anti  (24)  except  that  the  value  of  c(x),  required  to  determine 
d(x).  is  t he  fiat  plate  value  given  by  Eq.  (23)  with  the  curly  brackets  equal 
to  one  (since  0  is  equal  to  y.oro).  The  values  so  calculated  are  designated 
with  *.  As  before  we  let  the  prescribed  f/c  be  R  and  the  calculated 
f/c  be  U(n\  This  time  we  do  not  iterate  V*  and  p*  for  each  value  of 

•*  *  4 

T  (and  lienee  R  )  because  they  arc  loss  sensitive  than  R  to  changes 
in  T  .  It  is  also  necessary  to  adjust  the  chord  values  c^  (x)  and  (x) 
at  each  iteration  so  that  c  (x  =  1  )  =  l  -  R  and  c  (x  =  1 )  *  R.  This  adjust¬ 
ment  in  the  chord  values  is  necessary  because  the  6  values  are  defined 
only  in  the  regions  0  v.  c  >.  (1  -  R)  and  0<c  <_  R.  With  these  exceptions 
this  iteration  is  carried  out  in  the  same  way  as  in  the  flat  plate  iteration 
until 

v.(« -  1 )  .  v.(»> , .  ,„*(»- l).  ,*(«)  ,  i=d  ,  R*(a) .  R 
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are  less  than  the  allowable  error.  We  then  calculate  the  exact  value  of 
c(x}  (with  the  curly  brackets)  and  recalculate  V^»  and  R^.  If 

necessary  this  last  step  is  iterated  until 

j  v(n  "  1 )  _  y(n)  i  j  p(n  -  1 )  _  p(n)  j  aad  j  R(n)  _  R  j 

are  less  than  the  allowable  error.  When  the  camber  is  small  it  is  found 
that  it  is  not  necessary  to  iterate  this  last  step. 

Our  program  is  written  so  that  the  6  values  may  be  obtained 
either  by  a  four  point  interpolation  from  a  stored  set  of  6  versus  c 
values  or  from  an  analytic  expression. 


"* 1  a  s-e. 


Kig*  2  The  coordinate 


system  and  notations  for  a  flat  plate  with  a  flap. 


Fig.  3  The  variation  of  C.,  and  Cw  with  <J,  o  and  t  it  for  a  flapped  flat  plate  with  f/c  =  0.2. 


Continued  Fig.  7  Continued 


for  the  flapped  flat  plate  with  cr=  0,  The  solid  lines  are  lines  of  constant  a  and  the  dashed 
of  constant  <  K 
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The  pressure  distributions  for  a  flapped  flat  plate  with  f/ c  =  0.  2,  a  =  10  and  «  it  =  0  ,  20  ,  40  and  60 
compared  with  the  experimental  data  of  Meijer. 


Continued 
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Fig.  10  The  comparison  between  theory  and  experiment  for  and  C^.  The  open  symbols 
indicate  the  values  of  and  the  solid  symbols  those  of  obtained  by  integrating 
curves  through  the  experimental  data  in  Fig.  9. 
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Fig.  11  The  comparison  between  theory  and  experiment  of  for  a  flat  plate.  The  <7's  in 
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Parkin's  experimental  data  are  obtained  from  measured  cavity  pressures. 


comparison  between  theory  and  experiment  of  C-.  and  C_.  for  a  flapped 


Fig.  13  The  comparison  between  theory  and  experiment  of  C  r  and  for  a  flapped  flat  plate  with  f/c  =  0,  5  and 
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